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Abstract 

The mean value of the normalized current helicity a p — B ■ (V x B)/_B 2 in solar active 
regions is on the order of 10 _8 m _1 , negative in the northern hemisphere, positive 
in the southern hemisphere. Observations indicate that this helicity has a subsurface 
origin. Possible mechanisms leading to a twist of this amplitude in magnetic flux tubes 
include the solar dynamo, convective buffeting of rising flux tubes, and the accretion 
of weak external poloidal flux by a rising toroidal flux tube. After briefly reviewing 
the observational and theoretical constraints on the origin of helicity, we present a 
recently developed detailed model for poloidal flux accretion. 
Keywords: Sun: active regions, Sun: magnetic fields, MHD 



1 Introduction 

Three-dimensional vector fields may possess chirality (handedness), i.e. their 
structure may not be mirror symmetric. Concentrating on the magnetic field 
B, the simplest way to construct such a structure involves the superposition of a 
magnetic flux tube and a current tube. The result is a twisted tube where both 
the magnetic field lines and the current lines have a helical structure (Fig.l). 
This suggests that the helicity can be characterized by the scalar product of B 
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Figure 1: As Maxwell's equations are linear, the solutions can be superposed. The 
superposition of a straight flux tube (top left) and a straight current tube (top right) 
results in a helical flux tube (bottom). 



and j; as the latter is proportional to V x B, the current helicity is defined as 

h c = B • (V x B) (1) 

Another quantity characterizing twist, often preferred for theoretical calcula- 
tions is the magnetic helicity: 

h M = A • B = A • (V x A) (2) 

(To be precise, the above defined he and /im are helicity densities, helicity being 
their integral over a finite volume.) 

Twist is just one form of helicity. As illustrated in Fig. 2, the twist of a tube 
or band can be converted into writhe or kink: 




Figure 2: Twist and kink of a magnetic flux tube can be converted into each other, 
as easily demonstrated by pulling the ends of a kinked paper tape, then pushing them 
together again. 
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2 Observations of helical structure on the Sun 

The apparent morphology of solar atmospheric structures, espe cially on H„ 
ima ges, is often vorte x-like, twisted or helical. Early claims by lHald (|l927li 



and iRichardsonl ([19411 ) regarding the predominance of one chirality or another 
on a given hemisphere were hard to confirm owing to projection effects, a lack 
of information on 3D structure and large statistical scatter. In any case, the 
existence of a hemispheric rule of chirality was convincingly demonstrated at 
least in the case of quiescent prominences and interplanetary magnetic clouds 



( Zirker et alj 119971 ). Nevertheless, only direct magnetic measurements could 
provide convincing direct evidence for chirality rules in active regions. 

Most determinations of current helicity rely on a constant-a force-free mag- 
netic field model, constructed using the measured line-of-sight field component 
as lower boundary condition. Such models are based on the assumption that 
V x B = a p H where a p is constant. Multiplying this relation by B, it is clear 
that a p = he I B 2 is a normalized measure of current helicity, with a dimension 
of 1/length. Its meaning is easy to visualize, as l/a p is the length along a flux 
tube over which he field lines make a full circle around the tube. 

Of the possible values of a p , a best fit is then chosen by comparing the re- 
sulting magnetic field structure to either mor phological features seen in solar 



images taken in H a , EUV, etc. ( Seehafeilll990l ) or, more recent ly, to the hor- 



izont al field components as measured by vector magnetographs (jPevtsov et al.l 
19951) . A third method consists in simply taking the ratio tt p = (Vx T£) n /B n , as 



measured by vector magnetographs, and averaging it over th e active region. Re- 
sults derived by different m ethods generally agree quite well ( Leka fc Skumanichl 
ll999HBurnette et alJliool . 

Currently all these methods yield a single mean value of a p over the whole 
active region. This is because the low S/N ratio of vector magnetograph mea- 
surements inhibits a reliable study of current helicity distribution across the 
plage; for the same reason, usually only pixels with relatively high field strength 
(B > 100 G or so) are used in the determination. It is to be hoped that future 
improvements in magnetograph sensibility will remedy this situation. 

These studies have now firmly established that active region magnetic fields 
have helical structures, with a higher occurrence of negative helicity in the 
northern hemisphere. The observations show that the typical average value of 
the current helicity parameter a v in an active region is on the order of 10 -8 m _1 
van Driel-Gesztelvi et al. 20031) . 



The tilt of the axis of bipolar active regions relative to the E-W direction 
is generally interpreted as evidence for writhe in the emerging flux tubes giv- 
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ing rise to active regions. The origin of this writhe is well understood: it is 
the c onsequence of the Coriolis force acting on the downflows in the flux loop 
legs (|D'Silva fc Choudhuril fl993h . The sign of current helicity implied by the 
observed tilts is positive on the northern hemisphere, i.e. opposite to the twist 
measured in the solar photosphere; its amplitude, however, remains well below 
the a p values quoted above. 



3 The origin of helicity 

The basic question regarding the origin of the observed magnetic helicities is 
whether they are generated after the emergence of the flux loop or the flux 
emerges already in a helical form. Shearing of the photospheric footpoints of 
magnetic loops due to differential rotation or smaller scale granular and su 



pergranular motions can, in principle, generate a helicity of right sign (jDeVore 



120001 ); however, the amount of helicity generated in this way seems to be in- 



suffic ient to compensate for helicity losse s in CME's (Ivan Driel-Gesztelvi et al.1 
2003). In addition, in an important paper [Pevtsov et alJ l|2003l ) studied the time 
development of the helicity in young, emerging active regions and found a good 
correlation between emergence rate and helicity increase. In particular, further 
increase of helicity ceases once the emergence of the loop, as measured by the 
increase in footpoint separation and area growth, comes to a halt. This obser- 
vation seems to settle the issue in favor of a subsurface origin of the observed 
twist. 

There are quite a few subsurface mechanisms which may naturally introduce 
twist in the structure of emerging flux loo ps. Such proposed me chanisms include 
helicity generation by the solar dynamo ( Seehafer et al. 20031) and buffetin g of 



the rising flux tubes by helical turbulent motions ( Longcope et al. Il998h . A 
further possibility is the effect of Coriolis force on flows in rising flux loops 
(|Fan fc Go ng 2000). As we mentioned above, this process is responsible for the 



generation of positive writhe in active region flux tubes in the northern hemi- 
sphere, and thereby for the o bserved tilt o f active regions. As magnetic helicity 
is conserved in ideal MHD (|Berger|[l999h , the same process should then also 



give rise to a twist of opposite sense, so that the net helicity remains constant. 
The amount of helicity generated by this process, however, is too low to explain 
the observed values of a p . On the other hand, a strongly twisted flux tube can 
develop a writhe by means of the kink instability. This mechanism should lead 
to a positive correlation of twist and writhe, in contrast to the helicity con- 
servation argument above. Indeed, recent observations bv lLopez Fuentes et al 
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( 2003h and Holder et al. ( 2004 ) indicate that this mechanism may be important 
in the case of at least some active regions with unusually high tilt. 




Figure 3: During the rise of a toriodal flux tube (here shown in cross section, hatched) 
through the convective zone, field lines of the weak external poloidal field may get 
wrapped around it. A tier lUhoudhurl \200& ). 

A fur ther possibl e theor etical explanation of the observed helicity was pro- 
posed by IChoudhuri ( 2003f) . who suggested that the poloidal flux in the solar 
convecti on zone (SCZ) ge t s wra pped around a rising flux tube, as sketched in 
Fig. 3. IChoudhuri et al.l ( 20041 ) later showed th at this mechanism give s rise 
to helicity of the same order as what is observed. IChoudhuri et al.l ( 2004( ) also 
used their dynamo model to calculate the variation of helicity with latitude over 
a solar cycle and found that the latitudinal distribution of helicity from their 
theoretical model is in broad agreement with observational data. 

If the magnetic flux in the rising flux tube is nearly frozen, then we expect 
that the poloidal flux collected by it during its rise through the SCZ would be 
confined in a narrow sheath at its outer periphery. In order to produce a twist in 
the flux tube, the poloidal field needs to diffuse from the sheath into the tube by 
turbulent diffusion. However, turbulent diffusion is strongly suppressed by the 
magnetic field in the tube. This nonlinear diffusion proc ess was studied in an un- 
twisted flux tube bv lPetrovav fe Moreno-Insertisl (119971). The model was subse- 
quen tly successfully applie d for sunspot decay dPetrovav fc van Driel-Gesztelvil 
Il997l) . In a recent paper (jChatterjee et al.l l2006h we extended this model by 
including the poloidal component of the magnetic field (i.e. the field which gets 
wrapped around the flux tube) and we studied the evolution of the magnetic 
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Figure 4: P/ots of B' z , B'^ and a p as functions o/£ for a rising flux tube (case A). The 
different curves correspond to the profiles of these quantities at the following positions 
of the flux tube: 0.7Rq (thick solid), 0.85-R© (solid), 0.9Rq (dashed), 0.957?© (dotted), 
O.98i? (dash-dotted). 



field in the rising flux tube, as it keeps collecting more poloidal flux during its 
rise and as turbulent diffusion keeps acting on it. 

Consider a straight, cylindrical, horizontal magnetic flux tube rising through 
the solar convective zone. As all variables in this model depend only on the radial 
distance r from the tube axis and on time, we study the wrapping of the large- 
scale poloidal field around the flux tube by considering a radially symmetric 
accretion of azimuthal field by the flux tube. A further complication is the 
expansion of the flux tube during its rise, due to the decrease of the external 
pressure. This expansion is assumed to be self-similar, so that the Lagrangian 
radial coordinate £ is related to the Eulerian radius r by £ = F(t)r, where the 
expansion factor F(t) was taken from thin flux tube emergence models. In the 
Lagrangian frame, flux density is rescaled as B' z = B z /F 2 and B'^ = B^/F. 
With these notations and assumptions, the induction equation takes the form 

eB, '=^~U^) (3) 



dt £d£\ n ae 
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Figure 5: Same as Fig. 4 but the field inside the flux tube is not allowed to decrease 
below 3B e at any height. (Case B) 



at 



= F< 



(4) 



(Note that the advection term appears in the <f> component only, as it does 
not represent a truly radially symmetric inflow; instead, it is just designed 
to mimick the e ffect of the wrapping of poloidal fi eld lines around the rising 
tube.) Following IPetrovav fc Moreno-Insertid ( 1997I ) . the magnetic diffusivity r\ 
is specified in the form 

V = -, ; ^ 9 (5) 



1 



{B/B eq y 



All our calculations are performed for a flux tube of flux 10 22 Mx and initial 
field strength 10 5 G. 

Numerical solutions of the equations are presented in Figures 4 and 5. The 
conclusions drawn from our model hinge on some assumptions, especially con- 
cerning the subsurface magnetic field structure in the last phases of the rise of 
the tube. The field strength in the rising tube, as calculated from thin flux tube 
emergence models, decreases well below the turbulent equipartition value B e 
near the surface (i? 2 /2/i = pv 2 /2). The presence of 3000 G magnetic fields in 
sunspots is a compelling proof that magnetic fields may never fall to such low 
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values; in fact, at least in photospheric layers, they are at about 3B e . This is 
presumably the result of flux concentration processes such as turbulent pump- 
ing and convective collapse. Thus, it may be more realistic not to allow the 
magnetic field to fall below ZB e . Figures 4 and 5 present results without and 
with this constraint, respectively. 

Inspecting the lower panels of Figures 4 and 5 one finds that the typical 
value of a p in the internal parts of the flux tube is of order ~ 10 _8 m _1 at a 
depth of O.85i? in both cases. However, as the flux reaches the solar surface, in 
case A the component spreads out due to diffusion and its gradient becomes 
smaller, reducing a p by about one order of magnitude. Only if the magnetic 
field inside the flux tube remains stronger than the equipartition field (the case 
B represented in Figure 5) is the B$ component unable to diffuse inside so that 
its gradient remains strong and a p is of order ~ 10~ 8 m _1 even near the surface. 
This suggests that our case B may be closer to reality, i.e. during the rise of 
the flux tube from O.9i? to O.98i? effective flux concentration processes are 
at work, keeping the field strength at a value somewhat above the equipartition 
level. 

Note, however, that in the calculations presented here the amplitude and 
sign of the poloidal field was assumed not to depend on depth. For alternative 
assumptions, significantly different current helicities may result, so the above 
conclusion should be treated with proper reservation. Details of the radial 
dependence of the poloidal field strength may strongly depend on the dynamo 
model. 

A more robust feature of the current helicity distributions, present in all the 
lower panels of our plots, is the presence of a ring around the tube with a current 
helicity of the opposite sense. This is clearly the consequence of the fact that 
on the outer side of the accreted sheath the radial gradient of the azimuthal 
field, and thus the axial current, is negative. This is an inevitable corollary of 
the present mechanism of producing twist in active regions. A rather strong 
prediction of this model is, therefore, that a ring of reverse current helicity 
should be observed on the periphery of active regions, somewhere near the edge 
of the plage. 

4 Conclusion 

One rather strong prediction of our model is the existence of a ring of reverse 
current helicity on the periphery of active regions. On the other hand, the 
amplitude of the resulting twist (as measured by the mean current helicity in 
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the inner parts of the active region) depends sensitively on the assumed structure 
(diffuse vs. concentrated/intermittent) of the active region magnetic field right 
before its emergence, and on the assumed vertical profile of the poloidal field. 
Nevertheless, a mean twist comparable to the observations can result rather 
naturally in the model with the most plausible choice of assumptions (case B). 

It is likely that the accretion of poloidal fields during the rise of a flux tube 
is just one contribution to the development of twist. Its importance may also 
be reduced by 3D effects: considering the rise of a finite flux loop instead of an 
infinite horizontal tube, the possibility exists for the poloidal field to "open up" , 
giving way to the rising loop with less flux being wrapped around it. It is left 
for later multidimensional analyses of this problem to determine the importance 
of any such reduction. In any case, the results presented above indicate that 
the contribution of poloidal field accretion to the development of twist can be 
quite significant, and under favourable circumstances it can potentially account 
for most of the current helicity observed in active regions. 
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